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to the amplification of the time-dependent parameter, ui p hys{k / a) , which may be the modified 
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dispersion relation in trans-Planckian physics. This squeezed quantum state is a nonclassical state 

that has no counterpart in the classical theory. We have considered the nonclassical states such as 
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the gravitational wave perturbation when the mode leaves the horizon. 
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I. INTRODUCTION 



Gravitational perturbations such as density perturbation and gravitational wave would 
have been created from quantum fluctuations during the inflation period when the wave- 
length of the perturbation was much smaller than the horizon size jl| . And if the wavelength 
stretches out to the Hubble radius during the inflation, the amplitude would be frozen out. 
The power spectrum of the gravitational wave, if the initial state is an adiabatic vacuum 
state (Bunch-Davies vacuum) j^J, gives a scale invariant one with n s = 1 and ut = at the 
horizon crossing time, which is defined by [3( 

P s = A s k n °~\ P T = A T k nT , (1.1) 

where A s and At are the amplitudes of scalar and tensor perturbations. 

These gravitational wave and primordial density perturbations may have undergone a 
squeezing process due to a time- varying spacetime background jjj. If the Hamiltonian is 
given by a free Hamiltonian term (aa^ + a)a) and an interacting term (aa + a' a') in the 
Heisenberg picture where a and are annihilation and creation operators, respectively, the 
squeezing occurs due to the interaction with a background. Although the squeezing can 
result from the interaction with the background spacetime, it is also possible to squeeze the 
vacuum state through a parametric interaction of the time-dependent dispersion relation at 
the subhorizon scale. This fact may be tested in the trans-Planckian physics. 



lorizon size, when dated back to the early 
Though the Planck scale physics has 
jo treat the trans-Planckian 
considering non-vacuum 



The wavelength that corresponds to the present 
universe, might be smaller than the Planck scale | 
not been completely understood, several ways were suggested 
effects such as modifying dispersion relations, uj phys (k/a) 0, (| 

initial state QJEJ, and modifying the space-time commutation relation motivated from the 
string theory In Ref. p|, a squeezed vacuum state as an initial condition was used 
to calculate the consistency relations (r = At/A s ) which is the ratio of the power spectra 
between scalar and tensor perturbations in inflation. 

In this paper we shall consider in trans-Planckian physics whether the initial state may 
evolve to a squeezed quantum state due to the parametric interaction which comes from the 
modified dispersion relation. So the gravitational wave in the inflation period may exist, in 
general, in a squeezed quantum state, not the vacuum state. This squeezing process may 
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slightly modify the power spectrum of the gravitational wave perturbations for superhorizon 
scales by the order of -S where M is the cutoff scale. We shall consider several types of 
nonclasscial states such as squeezed, squeezed coherent, and squeezed thermal states. The 
coherent states are in the borderline between quantum and classical states. The expectation 
value of the fields with respect to coherent states follow the classical equations of motion. 
So the coherent states may be a convenient tool to deal with the classicality of the quantum 
fluctuations. Further, it may be not excluded that the universe would have been very hot 
from the beginning. This means that thermal effects would play an important role during 
the history of the early universe. Thus, thermal states may have some relevance in dealing 
with the gravitational perturbations. 

Nonclassical states that have no counterpart in classical theory are widely studied in 
quantum optics. Nonclassical states of light exhibit the features such as sub-Poissonian 
photon statistics and squeezing 11 1. Quasiprobability distribution functions 3], which 
play a similar, but not exactly the same, role of a classical probability, are used to quantita- 
tively measure the nonclassical behavior of states. A particularly useful distribution is the 
Glauber-Sudarshan P distribution that is defined in terms of the density operator jl^ 

p = J d 2 aP(a)\a)(a\, (1.2) 

where \a) is a coherent state and a is a complex variable. If the P distribution of a state 
is more singular than a delta function and is not a positive definite, the state is classified 
as a nonclassical state 14] . The regularity condition for the quasiprobability distribution 
function is used to test whether states are either classical or nonclassical. 

The organization of this paper is as follows. In Sec. II, we introduce the squeezed 
quantum states and quantize the gravitational wave perturbation in the inflation period. 
The squeeze parameter is found from the mode solution in the trans-Planckian theory. In 
Sec. Ill, we calculate the two-point correlation function in various nonclassical states and 
obtain the power spectrum of the gravitational wave perturbation for the superhorizon scale. 
We finally conclude and discuss in Sec. IV. 
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II. QUANTIZATION OF SCALAR FIELDS IN DE SITTER SPACE WITH 
SQUEEZED QUANTUM STATE 



The gravitational wave perturbation in the inflation period will be quantized using the 
two-mode squeezed state formalism. The metric of the spacetime in conformal time is given 
by 

ds 2 = a 2 (T])[-dr] 2 + (7^ + h^dx^x 6 ], (2.1) 

where hij is a traceless, transverse perturbation and is assumed to satisfy \hij\ <C 7^. In the 
momentum space, h%j is expanded as 

If d 3 k 

kt ^ X) = aJ (2^)^^ e ^ fc(r/)elk ' X ' (2 ' 2) 
where ey is a polarization tensor and //& is a complex function. Due to the reality condition 
of h^, one has = 

The perturbed action for the gravitational wave is 



5S = I d 4 x—[ti lJ W - F(V, a)hijh %3 ], (2.3) 

(2.4) 



2 L ' " 
dr]d 3 k 



VkV-k ~ -^(PkP-k + PkV-k) + ((7) - u2 ( k ) I PM-h 



where a prime denotes the derivative with respect to the conformal time 77, and F(V, a(r))) is 
a function of spatial derivative, V, and the conformal time, and becomes u(k) after Fourier 
transforming. The last term in the second line, a modified dispersion relation, takes into 
account the trans-Planckian effect. It represents a time-dependent non-linear property when 
the size of the perturbation mode is much smaller than the Planck scale, 

u 2 (k) = aW phys (k/a), (2.5) 

but recovers the linear relation, uj ~ k, when the mode scale is greater than the Planck scale 

as in the standard inflationary model. From Eq. (|2.4|) . one obtains the conjugate momentum 

dC k , a! 
K k = jri- = Vk, (2.6) 

and the Hamiltonian density 

'Hk = ^kl^'-k + K-kHk - £fc> 

a' 

= TlkK-k + M 2 (k)fJ,kfl-k H (TTfc/U-fc + TT-fc/ifc). (2.7) 

a 



The Hamiltonian density leads to the equation of motion 



/4 + (u\k) - °pj fi k = 0. (2.8) 
To quantize the fields we introduce two linear time- dependent operators [3 13 

a k (v) = -i (v>Z(v) ~ ~^(v)j fa + i(pt(jl)^k, 

a' 
a 



a j _ k (v) = i yp'kiv) - —<Pk(ji) ) fa - i<Pk(v)^k- (2.9) 
hrodinger o per a 

00,0,13 



Here we treated £i k and w k as Schrodin ger oper ators. These operators satisfy the quantum 
Liouville-von Neumann equation 



t-^a k (r ] ) + [a k (r ] ),H k (r ] )]=0, (2.10) 

and the same equation holds for a k . It should be remarked that the eigenstates of these 
operators are exact solutions of the time-dependent Schrodinger equation, when (p k is a 
complex solution to the equation of motion 

^+(u 2 (k)- a pj^ k = 0. (2.11) 

The commutation relation of the fields in the momentum space at equal times 

\jik,Tt k ,]=i6(k + V), (2.12) 

leads to the usual commutation relation 

[a k , al,] = S(k — k'), (2.13) 

when the Wronskian condition 

VkVt-vW k = i (2-14) 

is imposed. We may thus interpret the a k and time-dependent annihilation and 

creation operator. These operators do not, however, diagonalize the Hamiltonian ()2.7j) . but 
provide the exact quantum states of the time-dependent Schrodinger equation |l7|. 

There are two ways of squeezing quantum states. First, we choose some preferred complex 
solution to Eqs. (j2.11|) and (|2.14j) . which may be chosen based on, for instance, the minimum 
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particle creation postulate 19] or the minimum uncertainty ^(|. Then a general solution to 
Eqs. (j2.11|) and (|2.14j) is the superposition 

<Pks(v) = uo¥k(v) + v (pl(r)), (2.15) 

with 

| Mo | 2 -|^o| 2 = l. (2.16) 

Here the subscript s denotes a squeezed state. We may write uq and Vq as squeeze parameters 

Uq = coshr, vq = e~ 1 ^ sinh r. (2-17) 

Second, a given state may be squeezed due to a parametric interaction, here, time-dependent 
spacetime background. The preferred solution ip k at rj can be written as the superposition 
of tp k and (fl at a different conformal time tjq. Note that <f k (r]) an d w k {r]){= v'ijl) — ^^(v)) 
satisfy the Hamilton equations. Then the two functions 

Vk(v,Vo) = -i[&k(v)<Phto) — <Pk(v)™k{vo)], (2-18) 

satisfy Eq. (|2.11|) at rj and rj , and have the initial value u k (r] ,T] ) = 1 and v k (r) ,r) ) = 0. 
The inverse transformation of Eq. ()2.18|) is found to be 

¥ks(v) = w fc (77, 77o)^(%) +v k (v,Vo)<ft(Vo), 

r^ksiv) = u k (r], r]o)^k(r]o) +Vk{v>Vo)™k(Vo)- ( 2 -19) 

In fact, <fks(v) an d ^ k s(v) are solutions to Hamilton equations with the initial data ipk(Vo) 
and vj k (r) ). 

In the both cases, substituting Eqs. (|2.15|) and ()2.19|) into Eq. ()2.9j) . we obtain the 
Bogoliubov transformation 

aks(v) = u*a k (v) ~ v*a j _ k (ri), (2.20) 
al ks (v) = ual k (r)) -va k (rj). (2.21) 



In terms of the squeeze operator 

S k (z, rj) = exp [z k al(ri)al k (r]) - z* k d k (r])d^ k (r])] , (2.22) 
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where Zk = rke l< ^ k , the Bogoliubov transformation can also be written as 

ak S (v) = S k (z,rj)a k (r])Sl(z,T]), (2.23) 



-ks 



( V ) = S k (z, V )aUv)Sl(z, V )- (2-24) 



The squeeze operator S(z, rf) maps a number state into a squeezed number state which is 
an eigenstate of the squeezed number operator, Nks{v) = ^L( ? ?)^A:s(?7), 

\n,v)s = S(z,rj)\n,rj), (2.25) 

where 

\n,r]) s = \n,z;k,-k,i]). (2.26) 

We now solve the mode equation Eq. (J2.ll)) in each region during the inflation period. 
Consider an exponential inflation with the scale factor given by 

a( V ) = -± (2.27) 

When momenta k are greater than the cutoff scale, k c = aM, where M is the cutoff scale, 
the dispersion relation u(k) shows the non- linear "property. Several forms of uj n h,,Jk/a) in 

nnnna 

this region are suggested in literature |5J, |6J, UJ, la, |9J] . To solve the equation in this region, 
we assume that the adiabatic approximation holds, (|u/| u 2 ), and get the WKB-type 
solution 



If u phys (k/a) is given by 
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MV) = -^Le-J^W (2.28) 



uj phys {k/a) = Mtanh 1/p 



v 



^2.29) 



aM J 

where p is an arbitrary coefficient, the condition of the adiabatic approximation requires 
that | < 1. In fact, the small amplitude of cosmic microwave background fluctuations 
constrains -77 < 10 -5 while our present Hubble scale crossed the horizon during inflation 



2l\ . The mode crosses the cutoff scale when k = k c . If the wavelength of the mode is larger 
than the cutoff length scale, then the dispersion relation becomes linear in the physical 
wavenumber, u(k) — k, and leads to the exact solution 



The variances of the field and its conjugate momentum are calculated with the solutions 
(|2.28J1 and ()2.30|) to see whether the squeezing occurs during the expansion of the universe. 
In the trans-Planckian era, the variances with respect to the vacuum state are 

A/4 = (ji k n- k ) = -jj, (2.31) 
uj / 1 



An 2 k = (7r fe 7r_ fe ) = - ( 1 -I — ). (2.32) 

and 



2 y u) 2 rj 2 



1 

4 y ' uj 2 r] 2 



A^Att 2 = 7 ( 1 + _ . (2.33) 



Note that the uncertainty of the gravitational wave field hij and its conjugate momentum 
are related as 



A^=%i An^^Avr 2 , (2.34) 
a 



If to is given in the form of Eq. ()2.29|) . we can approximate to ~ k c = aM for k ^> k c . Then 
the uncertainty in the field, A/i fc , increases as the universe expands, whereas An k decreases. 
However, these have the minimum uncertainty as long as ^ < 1. So we can see that 
the squeezing occurs due to not only the expansion of the universe but the time-dependent 
parameter uj in the trans-Planckian era. On the other hand, for k < k c , the variances of the 
field and momentum with respect to the squeezed vacuum state are 

1 / 1 



= ^(i + dbr), (2-35) 



2k y ' k 2 r] 2 ' ' 
k 
2 



An% = * (2.36) 



and 



= \{ l + ]kf 2 )- (2 - 37) 

The modes whose wavelengths stay inside the horizon (krj 3> 1) have the minimum uncer- 
tainty. For the superhorizon scales (kr] 1), however, both the variance of the field, Afi k , 
and the uncertainty Afi k Airk increase with the conformal time. The scalar field in inflation 
can be described classically if the uncertainty Afi k Aiik is much larger than the minimum 
uncertainty [22]. So in this sense the gravitational wave can be treated as a classical when 
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the wavelength of the mode is much larger than the horizon length. The Bogoliubov coef- 
ficients u and v can be calculated from Eq. (j2.18j) where rjo denotes the conformal time in 
a region in which the trans-Planckian effect becomes important, and 77 is the time after the 
mode leaves the trans-Planckian cutoff scale. We choose the WKB solution (J2.28j) as the 
preferred complex solution (fikivo) i n Eq. (|2.19|) . which implies squeezing may occur through 
the trans-Planckian region. With the solutions ()2.28|) and ()2.30|) in each region, u and v 
take the form 



r / k 

^—ikrj+i J udrj I 



X 



1 



k 2 rf 



CO 

k 



% 

kr] 



1 + i 



\-±) J— L 

fa]) \krj kr) 0/ 



v(v,Vo) 



^—ikrj—i J iodr\ \ ^ 



x 



k 2 r] 2 k kr] ( k) ( kr]j \kr] kr]o / 



[2.38) 



[2.39) 



III. POWER SPECTRUM WITH NONCLASSICAL STATES 



In this section we calculate the power spectrum of gravitational wave perturbation when 
the perturbation mode stretches out to the Hubble radius during the inflation. The two-point 
correlation function at equal times 

dk sin kr 



(tt|/fa,x)/(77,x + r)|¥) 



k kr 



(3.1) 



gives the power spectrum defined by 

k a 



Pf(k) = — 2 J rf 3 re- ik - r (vl/|/(x)/(x + r)|vl/) 



2tt 2 
2^2 



(\h\ 2 )- 



(3.2) 
(3.3) 



Further, we consider the nonclassical states such as squeezed, squeezed coherent, and 
squeezed thermal states instead of the vacuum. Nonclassical states are quantum states 
that have no classical analog. These states have the quasiprobability distributions which 
may not always compatible with the interpretation as a probability distribution Q]. The 
quasiprobability distribution W(a, s) is used to show the nonclassical behavior quantita- 
tively and the existence of negative regions of distribution for nonclassical states. The 
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quasiprobability distribution, for example, the Glmmer-Sudarshaa P(o), is obtained from 
the complex Fourier transformation of Eq. (jl.2|) [12]. In Appendix El we briefly review the 
properties of the quasiprobability distribution, W(a,s). 



A. squeezed state 

First we consider the squeezed quantum state. In the previous section, we described the 
properties of the squeezed state. The ground squeezed state has a Gaussian wave function 



^o\^k,l^-k,v) = AAexp 



,w 



<Pk3 



(3.4) 



which is derived from the fact that a^s |0, rj) s — 0. Here M is a normalization factor. On the 
other hand, excited squeezed states (squeezed number state) have the non-Gaussian wave 
functions |23| 

M^) = (-1)" {^J Ln (j^p) ^\ (3-5) 



where Ln is a Laguerre polynomial. The temperature correlation functions of excited 
squeezed states show the nongaussianity. 

To calculate the power spectrum with the squeezed quantum state, we use the following 
relations 

s(a k a k >) s = u k vl(l + 2n k )5(k + k'), (3.6) 

s (al k al k ,) s = -u%v k (l + 2n k )6{k + k'), (3.7) 

s (a k a}_ k ,) s = \u k \ 2 (l + n k )5(k + k') + \v k \ 2 n k 5(k + k'), (3.8) 

s (a^_ k a kl ) s = \u k \ 2 n k 5(k + k') + \v k \ 2 (l + n k )5(k + k'), (3.9) 

where the inverse transformation of Eq. ()2.9j) was used 

a k = ua ks + v *aL ks , 

a]_ k = u*aJ_ ks + va ks . (3.10) 

Then the two point correlation function of hij with respect to the squeezed quantum state 
is 

>^IMx)^(x + r)K,> s = 1/ (2 f ) t /2(2 t)V2 ^ k, - (X+r) 
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PU k ) = TT^k 1 + 2n fe )[(K| 2 + \v k \ 2 ) + 2Ke(u* k v k )]. (3.13) 



x s (n,r]\iJ k (r))n k ,(r))\n,r]) s (3.11) 

x[|n fc | 2 + | Wfe | 2 + 2^eK^)]. (3.12) 

From this two-point correlation function, we can read off the power spectrum of the gravi- 
tational wave 

2n 2 a 2 

This gravitational wave power spectrum may be rewritten, using the squeezed mode (p ks in 
Eq. I|ZT5|) . as 

P 9W (k) = ^-J-^f(l + 2n k ). (3.14) 
For 7] > T) c , oj ~ k and kr/ c = ^, one has 

\u\ = cosh r ~ 1 H — — , 

1 1 4M 4 ' 

if 4 



v 1 2 = sinh 2 r ~ 



4M 4 ' 



1 i^ 2 f1M\ 

lZe(u*v) = - sinh 2r cos ~ cos f-^-j , (3.15) 

where we have used the mode crossing the cutoff scale at rj = rj c . The growing mode solution 
of Eq. (|2.11|) for k <C aH is given by (p k (rj) oc a{rf). Thus, the resulting power spectrum for 
the superhorizon sized perturbations is 



H 



P^)=[ — ) (l + 2n k ) 



2 



H 2 flM 



1 + — — r COS 



The power spectrum is corrected by the order |W and this result differs from Ref. |24j but is 



similar to Ref. 



M 2 \ H 



(3.16) 



in order of magnitude. This correction is too small to be detected by the 



CMB observations, but linear order corrections in (||) are expected to be detectable Q|. 
We use the result of Eq. (|A9J) to show that the squeezed vacuum state (n k = 0) with the 
squeeze parameter in Eq. ()3.15|) is a nonclassical state with the critical value (see Appendix 
A) 

Sc = e- 2 ^1-^. (3.17) 

This satisfies the condition for a state to be nonclassical (0 < s c < 1), which implies that 
the Glauber-Sudarshan P function (s = 1) is more singular than a delta function. 
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We may consider a more general initial state in the trans-Planckian period given by a 
squeezed state which is obtained from Eq. (j2.15|) 



(3.18) 



Here (p k (r] ) is the solution in Eq. (|2.28jl . Then the power spectrum for the superhorizon 
scales is 

H 2 ( ', , /2M\ 



P 9 M * (£) (l + 2n fc ) 



+ k>o cos 



2M 
IT 



2|u ||fo| cos 



2M 
IT 



(3.19) 



The terms containing the squeeze phase <p k may be set to zero due to the random character 
of the phase <p k of each mode (^tJ to yield 



H \ 2 

P gw {k) ~ ( —J (1 + 2n fe ) 



,2 , ,2 , ,2# 2 f 2M 

\uo\ + \v \ + k| Jp cos [-Jf 



(3.20) 



There is an increase of the overall factor by |uo| 2 + \vq\ 2 , which may be absorbed into 
the normalization of the spectrum. Once this normalization is done, the coefficient of the 
correction term can not be greater than one, thus slightly changing Eq. (|3. 16)1 . Note that 
the maximally squeezed states, r k — > oo, of the adiabatic vacuum have a factor 1/2. 

B. squeezed coherent state 



Now we consider the squeezed coherent state. The squeezed coherent state is defined as 
an eigenstate of the annihilation operators a ks {rj) and a_ ks {r)) 



a ks (r})\a,r}) s = a k \a,i]) s , 
d- ks (r])\a,r]) s = a- k \ac,r)) i 

where a k and a- k are complex constants, and 

|a,?7) 8 = \a k , a- k , 2,77). 



(3.21) 
(3.22) 

(3.23) 



The squeezed coherent state can be treated algebraically by introducing a displacement 
operator 



D(a) = D(a k )D(a^ k ), 



(3.24) 
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which is unitary 

t)(a)b\a) = D\a)b{a) = I. (3.25) 

The displacement operator translates by constant amounts the annihilation and creation 
operators through the unitary transformation 

b ] {a)a ks (r])b(a) = a k s(v) + a k, 

D\a)a_ ks (r))D(a) = a_ fes (r?) + a_ fc . (3.26) 

Similar relations hold for a ks and a_ ks . Hence one can see that the squeezed coherent state 
results from the unitary transformation of the squeezed state 

\a, V ) s = b(a)\ V ) s = b(a)S(z, V )\ri). (3-27) 

The expectation values of \i k and 7r fc with respect to the squeezed coherent state are 

fJ-kc = s (a,r]\{J,k\o£,r}) s = s {a,r)\(a k y k + a ] _ k (p*_ k )\a,rj) s 

= (u k a k + v*_ k a*__ k )ip k + (u*_ k a*_ k + v k a k )<p*_ k , (3.28) 
7i"fc c = s (a,rj\fj, k \a,rj) s = s (a, 77 1 I pf k - — fi k ) \a, r)) s 



a 



a' 



= /4 Wte- (3-29) 

a 

Note that Hkc(v) ^ s the solution of Eq. (|2.11|) . 

We calculate the following relations with respect to the squeezed coherent state 

s (a,rj\a k a k ,\a,ri) s = \u\a\ + v*_ k 2 a*_ k 2 ]5(k - k') + u k v* k 5(k + k') 

+2u k v*_ k a*_ k a k 5(k-k'), (3.30) 
s (a,T]\al k aL kl \a,r]) s = [u*_ k ai k + v 2 k a 2 k ]5(k - k') + u* k v k 5(k + k') 

+2u*_ k v k a*_ k a k 5(k - k'), (3.31) 
s (a,r)\a k aJ_ k ,\a,T]) s = \u k \ 2 5(k + k') + (u k u*_ k + v k v! k )a k a*_ k 5(k - k') 

+ [u k v k al + u*_ k v*_ k a*_ k 2 }5(k — k'), (3.32) 
s (a,rj\ai k a k >\a,r]) s = \v k \ 2 5(k + k') + (u*_ k u k + v * k v k )a*_ k a k 5(k - k') 

+ (u k v k a 2 k + u*_ k v*_ k a- k 2 )5(k - k'). (3.33) 
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With these relations, we find the two point correlation function of hij with respect to the 
squeezed coherent state 

1 r d 3 k 

s (a,77|%(x)/i^(x + r)|a,77) s = — J j^e^ 1 

x [M 2 {(K| 2 + \v k \ 2 ) + 2TZe(u* k v k )} + e 2 ^^ kc \. (3.34) 

We finally obtain the power spectrum 

Pgw = 2^[lv*r{(l«*| a + M*) + 2Ke(u* k v k )} + e 2ik ^ 2 c }. (3.35) 

The power spectrum consists of two parts; one is the contribution from the squeezed vacuum 
state as in Eq. (j3.13|) for n k = and the other is a classical part which is the expectation 
value with respect to the squeezed coherent state. The classical part \L kc satisfies the classical 
equation of motion f)2.8j) . As shown in Eq. (|A9|) . the regularity condition does not depend on 
the amplitude, a, of the squeezed coherent states so these states also satisfy this condition 
as in the case of squeezed vacuum states. 



C. squeezed thermal state 

We now use thermofield dynamics (TFD) for time-dependent system |28l l29j | and calculate 
the power spectrum of gravitational wave in a squeezed thermal state. The thermofield 
dynamics is a canonical formalism for finite temperature theory to describe quantum systems 
in thermal equilibrium. The idea of TFD is to double the system by adding a fictitious 
system and extend the thermal equilibrium in the system's Hilbert space to a thermal state, 
a pure state, in the extended Hilbert space of the total system. We briefly introduce the 
TFD formalism in one-mode state and after that calculate the power spectrum using the 
two-mode state TFD formalism. The density operators are defined by 

z 

p k (rf) = e -P"u&lmk(v) ^ (3.36) 

Zj 

where 

Z = Tie-! 3Hk . (3.37) 
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Here (3 and u are constants that may be fixed by the initial temperature and frequency. The 
density operator in the extended Hilbert space is given by 

Pk(v) = Pk(v) ® Pk(v) = ^e-^^W^W-aiW^W). (3.38) 

The thermal expectation value of the operator A of the system now takes the form 

(A) = Trp k (r])A = (0(/3), r]\A\0(P), 77), (3.39) 

where the thermal vacuum state is given by 

\0((3),ri) = y/l - e -^ e -(^AM&Uv)^ v ^ (3.40) 

with (0,77) = |0,6,77). As {ak,a\} and {a k ,a k } describe two independent systems, the 
nonzero commutation relations become 

[a k , a\,\ = 5(k - k'), [h, ~aU = 5(k - k'). (3.41) 

The thermal state is written as a time-dependent squeezed state of the vacuum state 

|0(/3),77)=f(^)|0,77), (3.42) 

where 

f{6) = exp[-0,(/3){a fe (77)a fe (77) - tfM"M}\- (3-43) 
Here 9 k (f3) is a temperature-dependent parameter determined by 

cosh0 fc (/3) = (1-e^)- 1 / 2 , 

sinh0 fe (/3) = e -^fc/2(! _ e -^ fc )-i/2_ (3 44) 

In the two-mode state formalism the thermal operator T{6) can be written as 

f(6) = cxp[-9 k (p){a ks a ks - a[ s a ks }} exp[-0 fc (/3){a_ fes a_ fes - al ks a_ ks }}. (3.45) 

Then the time- and temperature-dependent annihilation and creation operators through the 
Bogoliubov transformation become 

a k (P,r)) = f{6)a ks f\6) = cosh 9 k ((3) a ks (r)) - smh6 k ((3)~al k M, (3-46) 
~a k (P,v) = f{e)~a ks f\e) = cosh0 fc (/3)a fca (77) - sinh^(/3)a t _ fes (77). (3.47) 
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We get similar equations for a* k (j3, 77), a k (/3, 77) by using the Hermitian conjugate of these 
equations. The thermal state is a time- and temperature-dependent vacuum 

a k ((3, V )\0((3), V ) s = ~a k ((3,v)W),v)s = 0. (3.48) 

The thermal state |0(/3), 77) s , as an eigenstate of the invariant operators a k (f3, rj), is an exact 
eigenstate of the total system. To fix ui k in Eq. ()3.36j) . we assume that at rj — rji the universe 
was in thermal equilibrium state. And if we neglect the expansion of the spacetime at 77 = 77 j, 



the Hamiltonian is approximated as 

'Hk = K k n- k + u k fi k fi_ k , (3.49) 

to which we can set u k = a(rji)u p hys(k / a(rji)) . 

Using a similar method as in the previous sections, we get the following relations 

s (6(P),r)\a k a k ,\9(P),T]) s = u k vl(cosh 2 8 k + sinh 2 9 k )5(k + k'), (3.50) 

s (6(f3),7i\al k al kl \6(p),ri) s = u* k v k (cosh 2 6 k + sinh 2 6 k )5(k + k'), (3.51) 

s{e{(3),rj\a k al k ,\6{(3),r)} s = (|n fc | 2 cosh 2 9 k + |f fc | 2 sinh 2 9 k )S(k + k'), (3.52) 

s (9((3),r ] \al k a k ,\9(f3),r ] ) s = (\u k \ 2 sinh 2 9 k + |f fc | 2 cosh 2 9 k )5(k + k'). (3.53) 



Using the above relations, we can find the two-point correlation functions with the squeezed 
thermal state 

1 r d 3 k 

MP),v\h^)hH^ + r)\9(P), V )s = ~ 2 \ T^p^lrf 

x \{\u k \ 2 + \v k \ 2 ) + 2TZe(u k v k )](l + 2 sinh 2 fc )(3.54) 

from which we finally obtain the power spectrum 

P 9 M = ^2^K\ u k\ 2 + I^H + 2TZe(u k v k )](l + 2 sinh 2 9 k ). (3.55) 

Here 9 k follows from Eq. (|3.44jl as determined by the Bose-Einstein distribution at the initial 
time: 

s mn 2 ^ = -^±n- (3.56) 

Thus the resulting power spectrum with the squeezed thermal state is corrected by the ther- 
mal contribution factor (1 + 2 sinh 2 9 k ). The condition of the quasiprobability distribution 
W(a, s) being regular for the squeezed thermal state is from (|A9|) . 

•^TZT^^-wj- (3 - 57) 
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where we have used the relations ()3.15|) . In the limit of zero temperature, f3 — > oo, this 
condition reduces to that of the squeezed vacuum state. However, in the high temperature 
limit, (3 — > 0, this state may be regarded as a classical state. 

IV. CONCLUSION AND DISCUSSION 

We have studied the spectrum of gravitational wave perturbations in nonclassical states 
during inflation. Gravitational wave or primordial density perturbation, which is responsible 
for the large scale structure formation, would have been generated from quantum fluctuations 
when the wavelengths of the perturbation modes were much smaller than the horizon scales. 
After the wavelengths were stretched out to the horizon scale during inflation, they could 
be treated clS 8b classical perturbation. 

Any state can be, in the cosmological expansion, highly squeezed and this phenomenon 
is quite a quantum effect which cannot be explained in classical theory. Squeezing can occur 
when the time-dependent parameter is amplified at the subhorizon scale. We considered in 
this paper that an initial state can evolve to its highly squeezed state due to the amplifying 
time-dependent parameter which plays a role of the modified dispersion relation in trans- 
Planckian physics. So the gravitational field may exist, in general, in a squeezed quantum 
state during inflation. This is a dynamical process occurring through Eq. (j2.19j) . In Eq. 
(J2HEJ), the Bogoliubov coefficients u and v are functions of rj and rjo and satisfy the Eq. 




. From this point o 



Actually in Ref. 



view, the underlying physical concept differs from that of Ref. 



since the squeezed state was put as an initial condition, there 
was not a dynamical squeezing process. On the contrary, in this paper any initial state is 
squeezed through the dynamical process due to the time-dependent parameter. 

We also considered the nonclassical state whose behavior can be quantitatively described 
by the quasiprobability distribution such as squeezed quantum state, squeezed coherent 
state, and squeezed thermal state. The condition for the quasiprobability distribution to be 
regular is used to test whether the states are either classical or nonclassical. 

To compare with the observation, we calculated the power spectrum of the gravitational 
wave perturbation with the nonclassical states from the two-point correlation function. For 
the squeezed state, the power spectrum is corrected by (H) which is too small to be 



detected from the observation. It is expected for linear order, |j, corrections to be detectable 
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by the observations |26[. And the squeeze parameter satisfies the condition for the states 
being nonclassical. The power spectrum for the squeezed coherent state consists of squeezed 
vacuum state contribution part and classical part. The classical quantity /ifc c , which is the 
expectation value of the gravitational field with respect to the squeezed coherent state, 
satisfies the classical equation of motion (|2.8|) . The regular condition of quasiprobability 
distribution of the squeezed coherent state does not depend on the amplitude of the squeezed 
coherent state, so this state may be nonclassical as the squeezed vacuum state. And for 
the squeezed thermal state, the power spectrum is modified by the thermal contribution 
(1 + 2 sinh 2 8(f3k))- The nonclasscial condition depends on the parameter (3 and uok- 

Finally, the quasiprobability distribution measures the degree of the nonclassical behavior 
of states. This property will be used to investigate the transition from quantum to classical 
perturbations during the inflation for superhorizon scales. 
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APPENDIX A: QUASIPROBABILITY DISTRIBUTION 

We briefly review the properties of the quasiprobability distribution [l^ . The ensemble 
average of the operator F is represented by the density operator 

(F) = Tr(pF). (Al) 

This can be expressed in terms of the weight function W(a,s) by the integral form 

Tr(pF) = - / d 2 aW(a, s)f(a, -s), (A2) 

7T J 

where a is a complex variable. Though ^(a, s) can not exactly be interpreted as a probabil- 
ity distribution, it plays so closely the role of one that it is referred to as a quasiprobability 
distribution. The W(a, s) is identified with Glauber-Sudarshan (P), with the Wigner (W), 
and with Husimi (Q) distribution corresponding to the values s = 1, 0, and —1, respectively 
I^- 
IS 



Through complex Fourier transformations 3l|, W(a, s) is expressed in terms of the char- 
acteristic function s) as 



where 



W(a,s) = - fe^'- a '^s)d% 

7T J 



x(e,s) = Tr(pD(0)ef | « 1 



(A3) 



(A4) 



with the displacement operator D(£). The W(a,s) is normalized to unit through the trace 
of density operator 



and satisfies the Hermiticity 



1 = Trp = - / W(a, s)d 2 a, 

7T J 



W(a,s*) = W(a,s)*, 



(A5) 



(A6) 



which shows that the function s) is real for real values s. Especially, the P distribution 
is more useful in measuring the degree of the nonclassical behavior. If the P distribution 
of any state is more singular than a delta function and is not positive definite, this state 
is classified as a nonclassical state Q|. The P distribution no more singular than a delta 
function will be said to regular. 

The characteristic function of the squeezed coherent thermal states is derived 



X(f , s) = exp 



-A\£\ 2 - ^{B*£ 2 + BC 2 ) + C*i - CC 



(A7) 



where 



A 



1 -s 



+ n + (2n + 1) sinh 2 r, 



B = — (2n + l)e"^ sinh r cosh r, 
C = A, 



n 



i-l 



(A8) 



The condition for the integration of the quasiprobability distribution (jA3|) to be regular with 

0Q.Q 



the characteristic function (|A7j) is 



s < s c = (2n + l)e 



-2r 



(A9) 
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Then the quasiprobability distribution is calculated for squeezed coherent thermal states 



W{a,s) = (A 2 
x exp 



\B\ 



A 2 - \B\ 



A\a - Cf + \[B*{a - Cf + B{a* 



. (AlO) 



This result can be applied to the squeezed vacuum state by letting A = and n — 0, to the 
squeezed coherent state by n = 0, and to the squeezed thermal state by A = 0. Depending 
on whether a state satisfies the condition (|A9J) . it shows either a classical or nonclassical 
behavior. This result simply extends to the two-mode squeezed state if the characteristic 
function is given by the exponential of a quadratic form 



squeezed 

HQ- 
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